Effect of Quantum Fluctuations in an Ising System on Small- World Networks 
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We study quantum Ising spins placed on small-world networks. A simple model is considered 
in which the coupling between any given pair of spins is a nonzero constant if they are linked in 
the small-world network and zero otherwise. By applying a transverse magnetic field, we have 
investigated the effect of quantum fiuctuations. Our numerical analysis shows that the quantum 
fluctuations do not alter the universality class at the ferromagnetic phase transition, which is of the 
mean-field type. The transition temperature is reduced by the quantum fiuctuations and eventually 
vanishes at the critical transverse field Ac. With increasing rewiring probability, Ac is shown to be 
enhanced. 
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PACS numbers: 64.60.Cn, 46.60.Fr, 05.50.+q, 84.35. 

Various phenomena in nature and human society can 
be understood in terms of dynamical interaction among 
individual elements on a complex network. The equi- 
librium and dynamic properties of such systems depend 
strongly on the topology of the underlying network as 
well as how the individuals interact with each other. For 
example, the antiferromagnetic Ising model suffers from 
frustration effects on a triangular lattice, whereas it has 
a simple ground state on a square lattice. For the math- 
ematical simplicity, the network is often regarded as ei- 
ther completely random or purely regular. In reality, 
however, most of the biological systems, solid-state sys- 
tems, and human societies lie somewhere between these 
two extremes. Furthermore, it has recently been shown 
that such networks, now widely known as "small- world" 
networks, exhibit mixed properties: some of which are 
common to completely random networks, some common 
to purely regular networks, and others unique to small- 
world networks [1, 2]. 

A simple yet inspiring mathematical model for a small- 
world network has been proposed by Watts and Stro- 
gatz (WS) [3]. One starts from a regular network and 
"rewires" the links with probability p. As p varies from 
to 1, the resulting network "interpolates" from a purely 
regular network to a completely random network. The 
value of the model so generated comes from the fact that 
it captures important physics in a wide range of physical 
systems, which can be described by neither a regular nor 
a random network. 

In recent years, many authors have studied scal- 
ing properties, crossover behavior, percolation behavior 
[4, 5], the spread of infectious diseases [6, 7], signal- 
propagation speed [3], computational power, and syn- 
chronization [8, 9] of small-world networks. More re- 
cently, phase transitions of the Ising model [10, 11, 12] 
and XY model [13] on small- world networks have also 
been studied, where the crossover from one-dimensional 
to mean- field behavior has been found. Although most 
of the works mentioned above are concerned about clas- 
sical statistical problems on small-world networks, Zhu 
and Xiong [14, 15] have recently studied the localization- 
delocalization transition of electronic states on a small- 



world network which allows a quantum mechanical hop- 
ping among sites. 

In this work, we adopt a transverse-field Ising model 
on small- world networks, to investigate effects of the in- 
terplay between the unique topology of small-world net- 
works and the quantum fluctuations. From the previous 
works [10, 11, 12, 13], it is known that the small-world 
network topology enhances the correlation between the 
spins on the network and leads to the mean-field behav- 
ior of the system. Since the quantum fluctuations, intro- 
duced by the transverse field in our model, tend to de- 
stroy correlations, one can anticipate a non-trivial com- 
petition between the small-world network topology and 
the quantum fluctuations. 

Another motivation for studying the transverse-field 
Ising model on small-world networks is provided by the 
recent wide interest in quantum computing. A quan- 
tum computer can be regarded as controllable quantum 
spins (quantum mechanical two-state systems) interact- 
ing with each other through a network [16]. In realis- 
tic circumstances, the control of the spins, which can be 
achieved by carefully tuning the local magnetic field at 
each node and the coupling between each pair of spins, is 
imperfect. The effects of imperfections on quantum com- 
puting has been studied on a completely random network 
and have been shown to give rise to computation errors, 
which grow fast (exponentially or polynomially) with the 
number of quantum bits (or in short, qubits) [17]. It is 
therefore much worth studying the imperfection effects in 
terms of qubits on small- world networks. With the trans- 
verse field assumed uniform over the whole network, the 
model in our work may not directly represent a realistic 
quantum computer with imperfect controls, yet our work 
might provide a stimulation for studies in that direction. 

We consider N interacting spins in a uniform trans- 
verse magnetic field A. We assume that the interaction 
between spins is Ising-type. The Hamiltonian for the sys- 
tem is then given by 



H = H,+ Hr 



i<3 



(1) 
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where cr^, cr^, and cr^ arc the PauU matrices. The cou- 
pUng between the two spins at the zth and jih sites are 
defined on a small-world network (sec below for the pre- 
cise way the network is constructed) . Namely, we regard 
that the spins are placed on the small-world network, 
and that — J iii and j are linked on the network and 
Jij = otherwise [18] . The sign of J is irrelevant in this 
model since the ferromagnetic (J > 0) and the antifer- 
romagnctic ( J < 0) models arc interchangeable through 
a transformation that rotates every other spin about x- 
axis. Below, we will assume J > 0. 

The small- world network is constructed in a very simi- 
lar way as used by WS: We take a one-dimensional regu- 
lar network with 2k nearest neighbors [19] with a periodic 
boundary condition (o",|^jY = I = x,y,z). Among the 
total of Nk links we randomly choose pNk of them. One 
end of each chosen link is then rewired to a random site 
keeping the other end as a pivot point. Among the result- 
ing networks, wc discard those with disconnected parts, 
i.e., we only considers a network where all the nodes be- 
longs to a single cluster. 

Without going into details, several properties of this 
model may be qualitatively understood as follows. In 
the classical case (A = 0), this model is known to un- 
dergo a mean-field typo ferromagnetic phase transition 
at a finite temperature for arbitrarily small p > [10]. 
Finite transverse magnetic field introduces quantum fluc- 
tuations, which compete with the correlations enhanced 
by the topology of the small-world network. As a conse- 
quence, one expects A to suppress the transition temper- 
ature Tc- In the fully quantum mechanical case (T = 0) 
on the nearest-neighbor regular network (fc = 1), the 
Hamiltonian in Eq. (1) shows a quantum phase transi- 
tions at Ac = 1 of the Kosterliz-Thouless type [20] . 

Below wc will calculate the phase boundary, i.e., the 
transition temperature as a function of A, using the 
quantum Monte Carlo simulation. The universality class 
to which the phase transition belongs will also be deter- 
mined from various critical exponents. Before we give the 
results, we summarize the quantum Monte Carlo method 
specific to our model. 

We divide the inverse temperature (3 = ksT into M 
pieces with spacing e = (3/M. The partition function 
Z{P) = trexp {—[3H) is then approximated by the Trot- 
ter production formula 



Z(/3) « tr [exp {-eH^) exp {-eH,)] 



M 



(2) 



the error of which is on the order of {M oo). We 
now use the completeness relation 



(3) 



Si,t=±l 



where \Si^t) is the cigcnstatcs of the af, and put it be- 
tween the tth and {t+ l)th temperature shces in Eq. (2). 



0.8 



0.6 



0.4 



0.2 



□* 



X N = 200 
o 400 
+ 800 
□ 1600 
* 3200 



Ox 

" + >< 

*□+ o 



''*□□□+ O OOqo o 1 



2 

T/J 



FIG. 1: The fourth order Binder cumulant Un is drawn as 
a function of temperature for five different system sizes N. 
The parameters are given byfc = 2,p = 0.1, and B/J = 1. 
In the time direction, tlic size was fixed at Af = 30. As N 
grows, the crossing point converges to one single point, which 
is T/J = 2.04 in this graph. 



We may then write 



Z{0) = exp 

{S^.t} 



M 



where 



E E ^'^'■3 ^i'-t^o,t 

i<j t=i 

N M 
i=l f=l 

A = ^lncoth(eA) , 

B = ^ In [cosh(eA) sinh(eA)] . 



(4) 

(5) 
(6) 



Interpreting the temperature slices as imaginary time 
intervals, the above equation becomes a classical Ising 
model defined on a (1 -|- l)-dimensional lattice. The cou- 
pling in the spatial direction (indicated by the indexes i 
and j) is given by Jij and contains the topology of the 
small-world network. The coupling in the temporal di- 
rection (indicated by the index t), on the other hand, is 
a nearest neighbor coupling. 

Due to randomness in the topology of the network, 
this problem is naturally suited to a numerical analy- 
sis. We have thus performed a quantum Monte Carlo 
analysis using the above partition function. In order to 
increase efficiency near and below the ferromagnetic tran- 
sition temperature, a cluster algorithm was adopted [21]. 
Since our model is different from a usual regular Ising 
model, wc have developed a new cluster algorithm which 
takes into account the small-world network topology and 
the anisotropy between space and time axes. 

The transition temperature Tc(A) has been deter- 
mined from a finite size scaling method. Varying tem- 
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perature, we computed the fourth order Binder cumulant 

[22] 



Un{T) 



Km')] 
[(m2)2] 



(7) 



for several different values of N. This quantity varies 
from one to zero while temperature is swept from zero 
(maximum order) to infinity (maximum disorder). The 
two different shapes of brackets, (• • • ) and [•••], denote 
the thermal average and the average over rewiring config- 
urations, respectively. A typical result is shown in Fig. 1. 
For large N, there is a single crossing point Tg. It is note- 
worthy that the result is independent of the number of 
time slices M once it exceeds a finite lower limit Mc, indi- 
cating that the correlation length in the time direction 
is a finite fraction of /3. When T/J > 1, the value of Mc 
is typically ~ 30 and it increases at lower temperatures. 



2.5'^ I- o o 



p = 0.1 
0.05 



1- 



0.5- 



2 

A/J 



FIG. 2: Pliasc boundary of a transverse-field Ising model on 
a small-world network. The small-world network was con- 
structed following Watts and Strogatz starting from the one- 
dimensional regular network with k = 2 and rewiring the links 
with the probability p = 0.1 and 0.05. 

Figure 2 shows phase diagrams in T-A space for k = 
2 and two different rewiring probabilities, p = 0.1 and 
0.05. First of all, we find Tc > in the absence of the 
transverse field, which agrees with previous results [10, 
11, 12]. When a small A is turned on, Tc remains finite 
although it decreases with increasing A. At a fixed A, 
we find that Tc increases with p, just as in the no field 
case. 

Extrapolating the phase transition line to T = 0, one 
may obtain the quantum critical transverse field Ac. This 
is an extension of the quantum critical point of the reg- 
ular Ising model [20]. From Fig. 2, one can clearly see 
that Ac increases with p. This implies that the more the 
Ising system is rewired, the more resilient it is to quan- 
tum fluctuations. 

We now characterize the ferromagnetic phase transi- 
tion by determining the universality class to which it be- 
longs. In the pure classical case (A = 0), it is known 



that the phase transition is mean-field like [10, 11, 13]. 
To address the question whether this is still the case for 
finite values of the transverse field A, we investigated the 
scaling behavior of the phase transition described above. 
First of all, the Binder cumulant was fitted to a scaling 
function of the form 



Un{T) = u(^{T-Tc)N^/'' 



(8) 



Due to the infinite-range nature of the small-world net- 
work, the above exponent V describes the divergence of 

coherence number Nc instead of the correlation length ^ 
[13, 23]. More explicitly, we may write 



Nc 



DC 



IT-Tc 



(9) 



near the transition. 

The other critical exponents have also been obtained 
from various physical quantities by fitting them to scaling 
functions. For example, the specific heat per spin is fitted 
to 



c(T) = N'^^-'c ((T - Tc)N^^''^ . 



(10) 



Figure 3 shows an example of the scaling functions for 
(a) specific heat, (b) magnetization, and (c) susceptibil- 
ity per spin. If we choose appropriate exponents, results 
from systems of different sizes clearly collapse to one sin- 
gle curve near Tc- As T moves away from the scaling 
regime, the curves deviate. The best fitting values of the 
exponents are summarized in Table I. It turned out that 
they are the same as those of the mean-field transition 
to a very high precision. Therefore, we conclude that 
the quantum fiuctuations introduced by the transverse 
field do not alter the universality class of the ferromag- 
netic phase transition in the Ising models on small-world 
networks. 
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TABLE I: Critical exponents near the ferromagnetic phase 
transitions for A < Ac in an transver-field Ising model on a 
small-world network. For comparison, mean-field (MF) values 
are also provided. 

So far, we have studied the scaling properties of an 
Ising system on small-world networks at finite tempera- 
tures. However, the quantum critical point Ac at T = is 
in itself of much interest, since in general the universality 
class of a quantum critical point is different from that of 
classical transitions at finite temperatures [20] . Whether 
the small- world network will change the universality class 
of the quantum critical point as compared to that in a 
regular Ising model is a highly intriguing question. In 
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order to address that question, however, one has to use 
a different technique than those used above, because the 
results from the quantum Monte Carlo simulations be- 
come unreliable near T = 0. Therefore we leave it as a 
topic for further study. 

In summary, we have used quantum Monte Carlo simu- 
lations to obtain the phase diagram of an Ising system on 
small- world networks in the presence of a transverse mag- 
netic field. The ferromagnetic phase persisted at finite A, 
although the effect of quantum fluctuations introduced 
by the transverse field was manifested by the decrease of 
Tc. At a fixed A, we have also shown that Tc increases 
with the rewiring probability p. From various scaling ex- 
ponents, we have argued that the ferromagnetic phase 
transition at finite field was still of a mean-field type. 
Eventually, decreases to zero at a quantum critical 



point at a finite field Ac, but Ac increases with increas- 
ing p. Since the ferromagnetic region increases with p 
in both T and A directions, we conclude that the small- 
world topology of the spin system competes against both 
thermal and quantum fluctuations and enhances the cor- 
relation and ordering of the spins. 



Acknowledgments 

We thank M. Y. Choi, H. Hong, and B. J. Kim for very 
useful discussions. We acknowledge the support from the 
Swiss-Korean Outstanding Research Efforts Award pro- 
gram (SKORE-A). 



[1] A.-L. Barabasi, Linked: The New Science of Networks 

(Perseus Publishing, Now York, 2002). 
[2] M. Kochen, cd.. The Small World (Ablcx, Norwood, NJ, 

1989). 

[3] D. J. Watts and S. H. Strogatz, Nature 393, 440 (1998). 
[4] M. E. J. Newman and D. J. Watts, Phys. Rev. E 60, 

7332 (1999). 

[5] C. Moore and M. E. J. Newman, Phys. Rev. E 62, 7059 
(2000). 

[6] M. Kuperman and G. Abramson, Phys. Rev. Lett. 86, 
2909 (2001). 

[7] C. Moore and M. E. J. Newman, Phys. Rev. E 61, 5678 
(2000). 

[8] M. Barahona and L. M. Pecora, Phys. Rev. Lett. 89, 

054101 (2002). 

[9] J. Ito and K. Kaneko, Phys. Rev. Lett. 88, 028701 (2002). 
[10] A. Barrat and M. Weigt, Eur. Phys. J. B 13, 547 (2000). 
[11] M. Gitterman, J. Phys. A 33, 8373 (2000). 
[12] C. P. Hcrrcro, Phys. Rev. E 65, 66 110 (2002). 
[13] B. J. Kim, H. Hong, P. Holme, G. S. Jeon, P. Minnhagen, 

and M. Y. Choi, Phys. Rev. E 64, 56 135 (2001). 
[14] C.-P. Zhu and S.-J. Xiong, Phys. Rev. B 62, 14 780 
(2000). 



[15] C.-P. Zhu and S.-J. Xiong, Phys. Rev. B 63, 193405 
(2001). 

[16] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambride University Press, 
Cambridge, 2000). 

[17] P. H. Song and D. L. Shepelyansky, Phys. Rev. Lett. 86, 
2162 (2001). 

[18] A more realistic model on solids, Jij may also depend 
on distance between i and j. Here we adopt the simplest 
model. 

[19] We restrict A; > 2 since the WS model with fe = 1 has zero 

clustering coefficient and is exceptional in many ways. See 

for example, Ref. 10. 
[20] S. Sachdev, Quantum Phase Transitions (Cambridge 

Univ. Press, Cambridge, 1999). 
[21] M. E. J. Newman and G. T. Barkema, Monte Carlo 

Methods in Statistical Physics (Oxford University Press, 

Oxford, 2001). 

[22] K. Binder and D. Heermann, Monte Carlo Simulation in 
Statistical Physics (Springer- Ver lag, Berlin, 1988). 

[23] R. Botet, R. JuUien, and P. Pfeuty, Phys. Rev. Lett. 49, 
478 (1982). 



(a) 




(T-T^)N /J 



FIG. 3: Universal scaling functions atfc = 2,p = 0.1, 
B/J = 1, Tc/J = 2.04, and M = 30: (a) specific heat, 
(b) magnetization, and (c) susceptibility. Each quantity is 
measured per site. The legend in (a) is common to all three 
figures. 



